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1 Introduction 



In view of the persistent interest in the high energy emissions from astro- 
physical compact objects P, we present a continuation of the study on 
the electron-positron pairs production in the presence of strong Bcr = 
m'^c^/{he) ~ 4.4 x 10^^ gauss), time varying magnetic fields [3 13 II]- As 
we have pointed out in our previous papers, the physical situation we have 
in mind is the creation of e~ — pairs around astrophysical compact ob- 
jects like neutron stars or black holes that can produce such strong magnetic 
fields 13 • In this picture, the time variation of the magnetic field is a con- 
sequence of the rotation of the compact object or its gravitational collapse 
following a supernova explosion. Also, the creation of the electrons and of 
the positrons is only an intermediate step. In fact, our final goal is to give a 
possible theoretical interpretation of the so-called gamma-ray bursts (GRBs) 
that are thought to be originated as a consequence of a supernova explosion 
around such astrophysical compact objects [T]. In our formulation GRBs 
are pulses of photons produced by the electrons and the positrons created 
by the time dependent stellar magnetic field through their annihilation [Zj 
and/or as synchrotron radiation Pj. Nevertheless, we are aware that the real 
physical situation is more complicated than our theoretical model for various 
reasons. In fact, on the one hand we treat just one process of pair production 
while, actually, other processes (the production of e~ — e"*" pairs by photons, 
see the review jH]) are present and can give dominant contributions. On the 
other hand, we stress much more "microscopic" features of the pair (photon) 
production such as selection rules that allow or prevent some particular ele- 
mentary processes and so on by neglecting "macroscopic" aspects that would 
be very difficult to be treated analytically (the macroscopic structure of the 
compact object magnetic field or the presence of already created particles). 
These aspects are, instead, privileged in other papers where the microscopic 
details are not brought out (see e.g. [1211111 1121 CBl)- In these papers, inspired 
by the suggestion in (14j| that nonthermal radio emissions can be accounted 
by the gravitational collapse of a massive star, a similar interpretation is 
given about the origin of a GRB as a consequence of a supernova explosion. 
In turn, the supernova explosion is thought to be caused by a combination of 
gravitational and magnetic factors and much more emphasis is given to the 
order of magnitude of the energy released or to other macroscopic features 
of the physical situation. 

Instead, till now, we have neglected the effects of the gravitational field 
created by the neutron star or by the black hole by performing all the cal- 
culations in a fiat spacetime. Even if there are situations in which this can 
be safely done [13 HEI , it is interesting to study what happens if the effects 
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of the gravitational field are taken into account. Actually, since the early 
work by Hawking on black holes thermal emission ^Tj, it has been widely 
demonstrated from a theoretical point of view that a gravitational field can 
be responsible by itself of particle creation (see Refs. ^H] for particle creation 
in the presence of stellar gravitational fields and Refs. for cosmological 
particle creation). The formalism used to deal with this subject is that of 
quantum field theory in curved spacetimes where the gravitational field is 
treated classically and the consequent curvature of the spacetime metric is 
assigned, while the matter fields (Klein-Gordon or Dirac fields) are quan- 
tized starting from a Lagrangian density written in a general covariant way 
1^ 122] ■ In the same framework, the possibility of particle creation in the 
presence of gravitational and electromagnetic fields has also been investigated 

123 Ell- 
in the cited papers the particle creation is due to the fact that the gravita- 
tional field is time-dependent and the Bogoliubov transformation technique 
is used to calculate the production rate. Instead, as we will explain in the fol- 
lowing section, despite we start by using the same formalism, then we change 
our point of view. In fact, in our model the main process responsible for pair 
creation is still the magnetic field and the fact that it varies with time. The 
gravitational field is static and it enters the calculations because it modifies 
the one-particle states and energies of the electrons and positrons and con- 
sequently the production probabilities. In the present work, the structure of 
the gravitational field is assumed to be such that its effects can be calculated 
perturbatively. 

The paper is structured as follows. In section |21 the physical assump- 
tions and the consequent theoretical model is described: the gravitational 
effects are given starting from the Schwarzschild metric and then expanding 
it around a specific point, not too close to the event horizon and ending the 
expansion at the first derivative of the metric tensor. Then the magnetic field 
in this metric is given and the Hamiltonian of the second quantized electron 
field is built up. In section IHl the one-particle eigenstates and eigenenergies 
of the Dirac field in the given magnetic field with the perturbative effect 
of the gravity are displayed. In section |3] the pair production due to the 
time dependence of the magnetic field (in the constant gravitational field) 
are calculated with two explicit form of time variation. Some conclusions are 
presented in sectional Finally, four appendices contain explicit calculations 
that would have made heavy the main text. 

In what follows natural units {h = c = 1) are used throughout and the 
signature of the Minkowski spacetime is assumed to be H . 



3 



2 Theoretical framework 



As we have said in the Introduction, we want to calculate the probability of 
producing an electron-positron pair in the presence of a strong, time varying 
magnetic field and of a gravitational field created by a neutron star or by a 
black hole. Even if we assume that the spatial structure and the time evolu- 
tion of both fields are given, the problem formulated in these terms is very 
difficult to be faced. In fact, to determine a realistic form of the gravitational 
field and of the (electro) magnetic field we should fix the physical properties 
of the source (its mass, its electric charge, its angular momentum and so 
on) and solve the system formed by the Einstein equations and the general 
covariant Maxwell equations. Clearly, solving the Einstein equations coupled 
with the general covariant Maxwell equations is a hopeless problem and a 
number of approximations have to be done. In particular, we first assume 
that the Einstein equations and the general covariant Maxwell equations are 
disentangled. This corresponds to neglect the gravitational field produced 
by the magnetic field and to assume that the spacetime metric is determined 
only by the stellar object. But, in order to fit our requests this stellar object 
should be capable to produce a time-varying magnetic field. In this sense, 
it could be a charged rotating black hole or a magnetized collapsing neutron 
star but the corresponding spacetime metrics would be still very difficult to 
deal with. The problem can be definitely simplified by assuming that the 
spacetime metric is actually that produced by an uncharged, non rotating 
star and that the corrections to this metric due to the charge or to the rota- 
tion of the star itself can be neglected. In this approximations, our starting 
point is the metric tensor corresponding to the field created by a spherical 
body of mass M outside the body itself. If we call t, X, Y and Z the four 
coordinates, this metric tensor can be written as 



g^y{X,Y,Z) = diag 
where 



(ii.i) " 

F± (X, Y,Z) = 1± = (II.2) 

with Tg = 2GM the gravitational radius of the body {G is the gravitational 
constant). We have chosen the so-called isotropic metric instead of the usual 
(and equivalent) Schwarzschild metric, because from Eq. pi.l|) we see that 
the spatial distance is proportional to its Euclidean expression and this will 
simplify our future calculations. We point out that in this metric the event 
horizon of the body is the spherical surface Va^+T^ + Z^ = rg/4. 
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Now, an e~ — pair is created in a spatial volume with typical length 
given by the Compton length A = 1/m (m is the electron mass). Since A 
is much smaller than the gravitational radius Vg of a neutron star or of a 
black hole, we are allowed to make some simplifications on the metric tensor 
pLljl . In particular, we can always assume that the pair is created in a small 
neighborhood of the space point Pq labeled by the coordinates {Xq, 0, 0) with 
Xo > Tg/A + A and A > 0. If P = (X, Y, Z) = (Xq + x, y, z) is a generic 
point in this neighborhood then |a;| < A, \y\ < A and l^l < A and we can 
approximate the metric tensor g^v{P) calculated in P with the metric tensor 



dX 



X + 



dY 



y 



dZ 



P=Po P=Po P=Po 

{11.3) 

in which only the terms up to first order in x/Xq, y/X^ and z/Xq have been 
kept. It can easily be seen that gli^{x,y,z) actually depends only on x and 
that it can be written in the form 



g^Xx) = g'-^J + h,Ax) 



where 



(11.4) 



gjj^J = gfiuiPo) = diag 



dg,AP) 



dX 



X = diag{2gtx, 2gsX, 2gsX, 2gsx) 



P=Po 



with [see Eqs. (ITLI|1 and (jTL2 



(11.5) 
(11.6) 



4Xr 



1 + 



4Xn 



gt 



AXoJ 2Xo^ 



1 + 



4Xo 



1 + 



1 + 



4Xr 



4Xo/ 2XA 



:il.7) 

:ii.8) 



It is evident that, in order that gj}J{x) is a good approximation of g^i,{P), 
Xq can not be chosen to be too close to the critical value Vg/A. Just to give 
an idea, it easy to see that, if X ^ 1 is a large pure number, then 



(1), 



g^l^l{P) - gU 



X] 



g,,{P) 



< 



N 



no sum 



with /i = 0, . . . , 3, if 



Xn > 



XA. 



(11.9) 



(11.10) 
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This condition automatically implies that 



< — 1= no sum (11.11) 

2Vn 



(0) 

with /i = 0, . . . , 3 and then that h^y{x) can be considered as a small correction 
of g^l^} . In what follows, we assume that the previous inequalities hold with 
sufficiently large A^. 

Now, we pass to the description of the structure of the magnetic field 
that we will denote as B(r, t) and that, for simplicity, we will assume to lie 
in the y — z plane: B(r, t) = [0, -Bj^(r, t), i?^(r, t)]. With similar arguments 
used in the case of the gravitational field we can assume that the magnetic 
field is uniform and slowly varying in time in the spacetime volume where an 
e~ — e"*" pair is created. By writing it as 

B(r,t) = B(t) = I By{t) \=B{t)\ sin^9(t) I (11.12) 




COS'i9(t) 



with 



B{t) = ^lBl{t) + Bl(t), (11.13) 



tan^(t) = 1^ (11.14) 

and by also reminding that we work in the strong-field regime, the magnetic 
field is such that 

2 

TfL 

Bit) > Bcr = — , (11.15) 

e 

< ^ = m (11.16) 



B{t) X 

where — e < is the electron charge. Before going on we want to make two 
observations about the magnetic field. Firstly, we have shown that the mag- 
netic field pi.l2|l does not satisfy the Maxwell equations in vacuum and in 
the spacetime with the metric pi.4|l , but an electric current is needed to pro- 
duce it. Secondly, we point out that the assumptions that the gravitational 
field is static while B()f:) is time varying are not contradictory even if the two 
fields are produced by the same source. From a theoretical point of view, 
this circumstance happens, for example, for a spherical body which collapses 
keeping its spherical symmetry and without rotating. In this case, in fact. 



6 



the gravitational field of the body is static because of the Birkhoff theorem 
|2ni while the magnetic field is found to be time- dependent from energy- 
conservation considerations. Obviously, the previous conclusion is correct if 
the gravitational field generated by the magnetic energy can be neglected. 

Since the spacetime is curved and the metric tensor is not simply r^^,^, we 
must pay attention in defining the vector potential A^{r,t) that gives rise to 
B(t). To this end, we assume that the three-dimensional components of the 
magnetic field B(t) define the spatial-spatial components of the full covariant 
electromagnetic tensor F^,^(r,t) that is 

F,2{t) = -F23it)=0, (11.17) 
F^3it) = -F,,{t)=Byit), (11.18) 
F2i(t) = -Fi2(t)=5,(t) (11.19) 

while the mixed or the full contravariant components are built by using the 
metric tensor pi.4|l . Now, by definition 



F^,(r, t) = A,.^{r; t) - A^.,{r, t) = A,,^(r, t) - A^,,(r, t), (11.20) 

then by means of Eqs. pi.l7|l - pi.l9jl . we can choose the covariant vector 
A^{r,t) as 

Ao(r,t) = 0, (11.21) 
A(r,t) = l[rxB(t)L, (11.22) 

A2(r,t) = i[rxB(t)]^, (11.23) 

A3(r,t) = i[rxB(t)]^. (11.24) 

These relations define a gauge analogous to the so-called "symmetric gauge" 
jl] in the Minkowski spacetime. Finally, it is convenient to define the three- 
dimensional vector potential A(r,t) = [Ax{r,t), Ay{r,t), Az{r,t)] as 

A(r,t) = -l[rxB(t)] (11.25) 

where the minus sign has been inserted to have d x A(r, t) = B(t). 

As we have done in our previous papers |2 El HI, we will calculate the 
probability to create the e~ — e+ pair from the vacuum by applying the 
adiabatic perturbation theory up to first order in the time derivative B(i:) of 
the magnetic field j2Ij ■ To do this, we have to build up the second quantized 
Hamiltonian of a Dirac field \E'(r, t) in the presence of the slowly varying 
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magnetic field pi.l2|) and in the curved spacetime with the static metric 
tensor pi.4|) and to determine its instantaneous eigenstates and eigenenergies 
|27j . The Lagrangian density of this system is given by |20] : 



2 



(11.26) 



The presence of the gravitational field is represented in Eq. flll.26|) by the 
following quantities 

^7«(x)^det[^W(x)], (11.27) 



7«^(x) =7"e«^(a:), (11.28) 
i 

v 

where 7" are the usual Dirac matrices satisfying {7", 7^} = 27]"'^ (assumed 
in the Dirac representation), ea^^{x) is a tetrad corresponding to the metric 
(ITL4|1 [231201 and where a"^ = i[7",7^]/2. We remind that in Eq. (ITL261) 
the adjoint field ^{r,t) is also defined as ^{r,t) = \E'"'"(r, t)7'^ in a curved 
spacetime. 

Since g^ilu{x) has been split as in Eq. pi.4j) with the matrix h^j_y[x) much 
smaller than the matrix gf}, we are allowed to keep in Eq. pi.26|) only the 
first-order terms in h^^{x). To do this, we need to calculate the quantities 
pi.27|) - pi.29|) up to first order. Concerning the determinant g^-^\x), we have 
immediately 

(7«(x) ~ (7(°)[l + /i(x)] (11.30) 

where 

g^'^ ^ det(^W) = -Ml (11.31) 
h{x) ^ h/^ix) = 2 - ^) X. (11.32) 

Also, being the metric tensor gliJ{x) diagonal, we can choose a diagonal 
tetrad with 

= ^= ^ 1 (1 - . (11.33) 

ef"(i) = , ' ~ ^ (1 + no sum (11.34) 
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(note that while the Greek indices run from zero to three, the Latin ones 
run from one to three). By means of this tetrad it can be shown that the 
connections r^(x) are already first-order quantities equal to 



dx 



(11.35) 



where e^a^^ is the diagonal zero-order tetrad with 



,(0)0 



no sum. 



(0) 



[11.36) 
[11.37) 



Now, the procedure to calculate the Lagrangian density pi.26j) up to 
first order in h^jj{x) is identical to that used in the weak- field approximation 
|28^ EH] and we give only its final expression: 



'1 -qex) 



2v v^t 
(1 - gpx) 



[m-i\idom) - (z9o^)7°^] + 

^-i%di + eA(v,t)]^- 
-^idi-eAi{v,t)]-i'^ 



[1 — gMx)rn^'^ 



where we defined the three couplings^ 



1 + 







3rg 


Qe = 




2X1 






9t 


9p = 










4>t 






9t 


9 m = 


4>s 


4>t 



4Xn 



2Xi 



2X1 



1 - 



4Xo 
4Xo 



1 + 



1 + 



4Xn 



4Xn 



(11.38) 

(11.39) 
(11.40) 

fll.41) 



^We introduced three couplings for later convenience because, actually, only two of 
them are independent. 
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Note that the modifications induced by the metric tensor pi.4|) in the La- 
grangian density pi.38|) are hnear in and Qs (obviously) but nonhnear in 
(pt and (ps- 

The definition of the Hamiltonian is a controversial operation in a curved 
spacetime [SOI ISD IS2] • We shall adopt the same definition given in (HB] for a 
scalar field that is 

(^, (9i^, ^, di^, n«, flW, r, t) = flW (5o^) + (9o^')nW-^W(^, d^-^, ^, d^^, I 

(11.42) 

where 



n«(r,t) = 

n(i)(r,t) = 



(11.43) 
(11.44) 



are the first-order conjugated momenta to the fields \E'(r, t) and \E'(r,t) re- 
spectively. In our case, 



n«(r,t) 
nW(r,t) 



ri(l - Qex) 



7"*(r,t) 



(11.45) 
(11.46) 



and then the Hamiltonian density pi.42j) becomes 

(1 - gpx) 



%*[z9, + eA,(r,t)]^- 



-^'[2 9,-eA,(r,t)]Y^J + 
+ (1 - gMx)m^!m^ 

(11.47) 

and it results independent of n'^^^(r, t) and n*^^)(r,t). Coherently, the total 
Hamiltonian is defined as 



(11.48) 



We point out that this definition can be shown to be equivalent in our case to 
the alternative definition which uses the energy-momentum tensor (see 1^01 
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for a detailed discussion about the relation between these two definitions). 
Now, 

(1 - gpxm^>)f^ = d,[{l - gpx)^>^'% - - gpx)% (11.49) 

then the Hamiltonian density pi.47|) is equivalent apart from a derivative 
term to the asymmetric Hamiltonian density (we use the same symbol to 
indicate it) 



^W(*,(9i^,^,r,t) 



+ ^-i%di + eA(r,t)][(l - gpx)%\ + 
+ (1 - 5(Afx)m^'^|. 

(11.50) 

This Hamiltonian density has the advantage that it can be written in the 
form 



jrW(*,(9i^,^,r,t) 



[l- gEx)^^n^^\v,-id,t)m (11.51) 



where \E'^(r, t) = \E'(r, t)7° and where we introduced the one-particle first- 
order Hamiltonian 



[1 — gpx)oL ■ [—id + eA(r, t)] + 



1 — gEX [ 2 

+ a ■ [-id + eA(r, t)]{l - gpx) + (1 - gMx)l3m 

1 



(1 - gpx)a ■ [-id + eA(r, t)\ + 



+ OL ■ [-id + eA(r, t)] (1 - gpx) + (1 - gMx)(3m+ 
1 



with /? = 7° and cx 



+ 9ex 



-.ct ■ [-id + eA(r, t)] + jSm 



[11.52) 



{P'y^ , P'y'^ , P'y^) . Before explaining the 



reason of this apparently unusual definition we observe that the one-particle 
Hamiltonian pi.52j) can be written as the sum 



n^^^ (r, -id, t) = (r, -id, t) + J(r, -id, t) 



(11.53) 



11 



of the zero-order Hamiltonian 



n^^\v,-id,t) = la - [-id + eA{r,t)] + V^/3m| (11.54) 

\/ (Da V J 



and of the first-order interaction 

9p 



J(r, —id, t) 



XOL ■ [—id + eA(r, t)] + a ■ [—id + eA(r, t)]x 
1 



- gMprnx+ 



+ 9eX 



a ■ [-id + eA(r,t)] + f3m 



^t{gp - gM)Pmx - ^ (r, -id,t)+ 7^ (r, -id,t)x] + 



+ gExn^°\r,-td,t). 



(11.55) 



In order to understand the definition pi.52|l we have to introduce the 
scalar product for spinors in our curved spacetime. If (^i{r,t) and (^2(1", ^) are 
two spinors it is defined as 



(Ci,C2)= / rfS^^-(7W(x)Ci(r,t)7«^(x)C2(r,t) 



(11.56) 



where S' is a Cauchy hyper- surface. We can choose S as the t = const, 
hyper-surface, then 

dSf,= {dr, 0,0,0) (11.57) 
and the scalar product ()II.56|) becomes 



(Ci,C2)= / dv^-gW{x)Cl{r,t)ei^^\x)C2ir,t) 



(11.58 ) 

where we used the definition flll.33|l . Finally, by exploiting Eqs. pi.3fljl - 
pi.32|) . the volume element can be written up to first order as 



and we obtain 



(Ci,C2) = I rfrV^(l-^?^x)Cl(r,t)C2(r,t). 



63(1 - gEx) 

(11.59) 

(11.60) 
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Now, with this definition of the scalar product we reahze that the "asymmet- 
ric" Hamiltonian pi.52|) is actually an Hermitian operator up to first-order 
terms. The fact that this definition of one-particle Hamiltonian is well posed 
is also corroborated by the form of the equation of motion of the field \E'(r, t). 
In fact, starting from the Lagrangian density (HLHHI), it can be seen that the 
equation of motion 

^ ^--^ +d = (11.61) 



dt d{do^) d{d^) 



can be written as 



ido"^ = n^^\r,-id,t)'^. (11.62) 
Finally, by using Eq. pi.Blll . the total Hamiltonian pi.48|l becomes 

H^^\t)= j dr^/^^{l- gEx)¥{r,t)n^^\r,-id,t)'^{r,t). (11.63) 

This Hamiltonian depends explicitly on time only through the dependence 
of T-l^^^r, —id,t) on the magnetic field B(t), then because of the condition 
pi.l6|) it is a slowly varying quantity. In the next section we will determine 
its instantaneous eigenstates and eigenvalues that will be used to compute 
the pair creation probability by means of the adiabatic perturbation theory 

123 

3 Determination of the instantaneous eigen- 
states and eigenenergies of the Hamilto- 
nian H^^\t) 

The first step to determine the instantaneous eigenstates and eigenenergies 
of the Hamiltonian pi.63|) is to quantize and diagonalize the corresponding 
constant Hamiltonian. If we introduce the constant magnetic field 



flH.ll 





( ° ^ 




f 




By , 




sin 











with 



B = ^IBI + Bl (111.2) 



tant9 = -^ (IH.3) 
By 
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and the vector potential in the symmetric gauge 



A(r) = --(rxB) 



(111.4) 



then the Lagrangian density ()II.38|) . the one-particle Hamiltonian pi.53|) and 
the total Hamiltonian pi.63|) become respectively 



+ 



(1 - gpx) 



+ eAi{r)]'^- 
-^^,-eA{r)]Y'^ 



- {1- gMx)m^'^j, 

(III.5) 

n^^\r, -id) = H^^\r, -id) + J(r, -id), (III.6) 
ijW = J dr^/^^{l- gEx)¥{r,t)n^^\r,-id)'^ir,t) 



(111.7) 



where [see Eqs. pi.54|) and pi.55|) 



n^'^\r, -id) = |a ■ [-id + eA(r)] + ^/^sPm\ 



(III.8) 



X(r, -id) = y^tigp - 9M)pmx - y [xH^^^r, -id) + -id)x] + 



+ gExn^'^\r,-td). 



(111.9) 



Also, the equation of motion pi.62j) becomes 

ido^ = n^^\r,-id)<I/ (III.IO) 

then, if 

^(r, t) = ^±,j{r) exp{TiWjt) (III.ll) 

with J a set of quantum numbers and \l/-|-j(r) are a complete set of modes 
with energies > 0, then 

H^^^r, -za)^±,j = ±Wj^±,j. (III.12) 
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After applying the charge conjugation operator to the positive-energy solu- 
tion \&_|_ ,,(r) of the equation piI.12|) with e instead of — e in the one-particle 
Hamiltonian T-l^^\r, —id) and calling the corresponding spinor Vj{r) (see ji] 
for details), Eq. piI.12j) splits into two equations: 

n^^\r, -id)Uj = WjUj, (III.13) 
n^^\Y, -id)Vj = -WjVj (III.14) 

where Uj{r) = \E'+j(r) and where we called Wj > the energies of the positron 
states that, in the presence of a magnetic field, have, in general, a different 
dependence on the quantum numbers j The states Uj{r) and V^(r) are 
assumed to be an orthonormal basis with respect to the scalar product pi.60|) . 
i.e. 

{U„Uy) = {V„Vy)=S,,y, (III.15) 
{U,,Vy) = iV,,Uy)=0. (III.16) 

Equations piI.13j) and piI.14j) with the orthonormalization relations piI.15jl 
and piI.16j) will be solved perturbatively in the next paragraph. At the mo- 
ment, we observe that, in order to quantize and diagonalize the Hamiltonian 
pil.7|l . we have to expand the Dirac field \l/(r, t) in the basis [Uj{r), Vj{r)] as 

^{r,t) = J2 [cAtmr) + dl{t)V,{r)] (III.17) 

with Cj(t) = Cjexp{—iwjt) and (ij(t) = (i]exp(iwjt) and we have to impose 
the usual anti-commutation rules among the coefficients Cj, cj, dj a nd (i| t hat 
are now operators. In this way, by using the equation of motion ()III.10|) . it 
can easily be shown that the Hamiltonian ()III.7|) assumes the usual diagonal 
form 

(in.18) 

where Nj = c^^Cj, Nj = rf^rfj and where we neglected the vacuum energy. ^ 
Obviously, if |0) is the vacuum state the eigenstates of this Hamiltonian are 
the Fock states 

^Note that to neglect the vacuum energy can be unsafe when we want to determine 
the time evolution of the gravitational field in the presence of quantum matter fields as 
its sources [50] . But in our case this fact does not cause any problems because the time 
evolution of the gravitational field is given (actually, in the present case the gravitational 
field is static). 
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with the corresponding eigenvalues 



E = J2^VJl+^5l^3l)■ 



In the next paragraph, we proceed by exphcitly computing the one- 
particle spinors Uj{r) and Vj{r) and the energies and and we will come 
back to the problem of the instantaneous eigenstates and eigenvalues of the 
Hamiltonian ()II.63|) in the paragraph ()3.3|) . By looking at the structure of 
the one-particle Hamiltonian pil.6|) . we see that the eigenvalue equations 
piI.13|) and piI.14|) can be solved perturbatively with respect to the inter- 
action Hamiltonian X(r, —id) by using the time-independent perturbation 
theory jHlj. In particular, we will compute the eigenstates Uj{r) and V^(r) 
and the eigenenergies Wj and Wj up to first order in the gravitational couplings 
Qe, gp and qm- 

3.1 Computation of the one-particle states Uj{r) and 
V^(r) up to zero order and of the energies Wj and 
Wj up to first order 

We first want to determine the zero-order solution of Eqs. piI.13|) and 
piI.14|l . Up to this order those equations can be written as 



These equations are the eigenvalue equations of a particle with charge — e 
and mass in the presence of the magnetic field ()III.1|) in the Minkowski 

spacetime, then we can write immediately their solutions We have 

already studied these solutions in |1] and we refer the reader to that paper 
for a detailed discussion about the meaning of the quantum numbers that 
will be introduced, their ranges of variation and so on.'^ Here, we only quote 
the results that will be useful later. Firstly, we remind that a complete set of 
commuting observables of this physical system is made by the Hamiltonian 
7^(0) (r, —id) , the longitudinal linear momentum V\\ , the longitudinal total an- 
gular momentum jTji = C\\ +S\\ and the transverse square distance R\ where 

■^We warn the reader that some unavoidable changes of notation have been done with 
respect to Ref. 0]. 
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"longitudinal" and "transverse" are considered with respect to the direction 
of the magnetic field. For this reason, the eigenstates and the eigenenergies 
are characterized by four quantum numbers that are j = {nd,k,a,ng}. Ac- 
tually, the eigenenergies of the electron do not depend on Ug and those of 
the positron do not depend on Ud [see Eqs. (15) and (34) in Ref. This 
circumstance leads us to embody in r the quantum numbers {n^, k, a} and 
in q the quantum numbers {Ug, k, a}. In conclusion, we have 

j = {nd,k,a,ng}, (III. 23) 

r = {nd,k,a}, (III.24) 

q = {ng,k,a} (III.25) 

and, consequently 

j = {r,ng} = {q,nd}. (III.26) 

With these definitions the electron and positron eigenenergies wf'^ and 
are given by the modified Landau levels 

= ' l — ^(j)sm^ + P + eB{2nd + 1 + a) = J <f)tm^ + 7^ + eB{2nd + l + 



(III.27) 




m2 + A;2 + eB{2ng + 1 - a) = J (ptm^ + ? + (^B{2ng + l-(x)]. 

(III.28) 

The corresponding eigenstates will be indicated as Uj{r) and Vj{r) respec- 
tively and, since we have already studied them in we will quote some of 
their properties in appendix A. 

As it is evident from the expressions ()III.27|) and ()III.28|) . the one-particle 
eigenenergies have two kinds of degenerations. The first one is due, as we have 
said, to the fact that the electron eigenenergies do not depend on the quantum 
number Ug and, symmetrically, the positron eigenenergies do not depend on 
the quantum number Ud. The second one is due to the fact that the electron 
eigenstates with quantum numbers {r+, Ug} = {rid, k, +1, Ug} and {r_, Ug} = 
{ud + l, k, —1, Ug} have the same energy whatever Ug and, symmetrically, the 
positron eigenstates with quantum numbers {q+,nd} = {nd, k, +l,ng + 1} 
and {q-,nd} = {ud, k, —l,ng} have the same energy whatever n^. This 
means, following the time-independent perturbation theory [34J, that the 
eigenstates Mj(r) and Vj{r) will not represent, in general, the correct zero- 
order eigenfunctions Uj^\r) and v}^\r). For this reason we indicated them 
by means of the symbols Uj{r) and Vj{r). 
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Now, we will compute explicitly only the zero-order electron eigenstates 
and the first-order electron eigenenergies, while the analogous results for the 
positron eigenstates and eigenenergies will be only quoted. Following the 
time-independent perturbation theory for degenerate states, we write the 
zero-order solutions of Eq. piI.13jl with a given energy as linear combinations 
of all the degenerate eigenstates Uj{r) and Vj{r) with that energy. We will 
characterize the new states by means of the index xq and this choice will be 
understood later: 



JjiO) ( ) ^ p{0) / N 



[Pr_'*,a;o;r_,ng'"r_,ng(r) + Pi},xo;r+,ng'^r+,ng{T') 

Ug 

^il,xo;r-,ng'^r-,ng{^) + ^il,XQ;r+,ng'^r+,ng{^) 



(III.29) 
(III.30) 
(III.31) 



In the first of these equations we considered that the eigenenergies of the 
so-called electron transverse ground states characterized by the quantum 
numbers {rQ^Hg} with tq = {0,k, —1}, are degenerate only with respect to 
Ug. Note that in Eqs. pil.30|) and pil.31|) the resulting spinors have been 
characterized, with an abuse of notation, by the quantum numbers r_ and r_|_ 
respectively, even if, in general, they will be only eigenstate of the linear lon- 



gitudinal momentum V\\ with eigenvalue k. But, the coefficients P^o,xo;ro,ng, 

m 

(III.32) 



Pr±,xo]r±,ng and P'"":^,xo;r±,ng, are the solutions of the secular equations 



>(°) and P(°) 

'±,xo]r±,ng ciiiu r r^,: 

{-^ro,ng;ro,n'g — ^ro,xo^ng,n'g) ^i(!,xo;ro,n' ~ 0' 



^r-,xo^ng,n'g^ P 



(0) 



r_ ,xo;r— ,n; 



r-,na\r+,n 



(0) 



g,' + ,"'g r_,xo;r+,7i' 



A 1 

r+,ng;r+,n'g tr_ ^xo^ng ,n'g J r-,xo;r+,n'g 



i""-g,n' ) ^r-,xo:r+,n' 



(0) 



|^(2^r_,?ig;r_,n^ ^r+,xo^ng,n'g) ^l^,: 



'r+,ng\r-,ng r-,xo;r-,n' 



■J p(0) 

-^r-- ,"3 ;»■+ ,"9 r+ ,xo ;r+ 



(-^'•+,ng;r+,n^ er+,i'0^ng,n^) Pr+,xo;r+,n' + ,'19;''- ,"9 P 



(0) 

r+,a:o;r_,n^ 



= 0, 

(III.33) 
= 0, 

(III.34) 

= 0, 

(III.35) 
= 

(III.36) 
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where, in general, 

= J dv^/^^u]{r)I{r,-id)uj,{r) (111.37) 

and where e^o^xo! ^r±,xo ^^^^ the first-order corrections (to be determined) to 
the eigenenergy of the eigenstates Uro]xo (r) and f/r± .xo (r) respectively. In ap- 
pendix B we show that the matrix elements 2r±,ng;r^,n'g vanish.^ This means 
that the perturbation X(r, —id) does not remove the energy degeneracy of 
the eigenstates characterized by {r_,ng} and {r+,np} and that all the co- 
efficients Pr±,a;o;r^,ng Can be put cqual to zero. In other words, we have to 
diagonalize the perturbation X(r, —id) inside every subspace labeled by the 
quantum numbers r = {nd,k,a}, and the zero-order eigenstates are indeed 
characterized by the quantum numbers r. Also, equations pil.29|) - pil.31|) 
become simply 

U^li^) = Y.^i%;r,ngUr,ng{r) (III.38) 

while equations ()l I l.32;i - ()l I l.3fi;i can be written as the single equation 

Now, by means of the same technique used in appendix B we can show that 
the matrix elements Ima-m' can be written as 



T 

-^r,ng; 



dxdye'* (x, y)xoe' , (x, y) 



{gp - 9m) + {qe - gp)w[^'^ 

Wr 

(III.40) 

where the operator xq and the functions 0'^^^^ {x,y) have been defined in 
Eqs. ffXT8|l and |X29|l respectively. Now, from^Eqs. (jXH), (jMI) and |X26|) 
we see that only the transverse functions d'n^ng{^->y) depend on 

Hg, then we have to determine the coefficients Pf}xo;r,ng in such a way that 
the linear combination Pr-°2o;r,ng6',^^ „^(a;, y) diagonalizes the operator xq.^ 
This linear combination is given in ^] [it is Eq. (104) in Complement Eyi]. 
The coefficients Pf}xo\r,ng result independent of the quantum numbers r and 
are given by 

PfL.ng = = \l^^^^ng{V^-0) exp (-'-^) (III.41) 



^Note that the matrix elements .„g -r^ and X,-^ .n' ;r_. .rig are not complex conjugated 
numbers because the operator X(r, —id) is not Hermitian with respect to the scalar product 
implicitly used in the definition ljIII.37|) . 

^Now, it is clear why we called the additional index xq. 
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where 



B.nAVeBxo) = ^ ds(V^xo + is)'''exp(-s^) (III.42) 

is the rigth-order Hermite polynomial In this way the spinors ?7i%(r) 
have the same form of the spinors Mj(r), but with the function 6*^^ „^(a;, 
substituted by the function 



Ug 

By using the expression ()A32|) of 6*^^^ we have 



(III.43) 



TT V 2 



-— I a^o + — ^ — 



e5 



■X 



TT 



(isexp (— s^) — r (VeBxo + is)VeB(x — iy) 



'nd\ 



TT V 2 



-— I % + — ^ — 



eB 
~2 



■X 



TT 



X exp [eBxo{x — iy)] / ds exp — s + eB{x — iy)s 



1 



V vr V 27r 



exp [(a; - xo)^ - iy{x - 2xo)] 

(III.44) 

where we used the formula p?7j 

/oo / 1.2 \ 

ds exp (-6?s2 ± = ^ exp ( Re(6?) > 0. (III.45) 

From now on, the quantum number Ug completely disappears in our cal- 
culations because it is essentially substituted by the quantum number xq. 
This circumstance allows us to simplify the notation. In fact, we can elimi- 
nate the subscript d from the remaining quantum number Ud and from the 
related operators, that is [see Eqs. ()A10|) and ()A11|) ] 



ad 




2 

eB 



{x - iy) + 





X + iy) 



d_ ld_ 

eB \dx i dy 

d_ _ ld_ 

eB \dx i dy 




Ud 



n. 



(III.46) 

(III.47) 
(111.48) 
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Also, later calculations will be simplified if we discretize the eigenvalues Xq. 
This can be done by imposing in the functions Q'n^^xoi^^y) periodicity 
condition at x = 

eBxoY\ f eBxoY 
— = exp —I 



exp [^—^J = exp [-^^^) (III-49) 

where Y is the length of the quantization volume in the ^/-direction. In this 
way, the allowed eigenvalues are given by 

2£7r 

xo;^ = — - i = 0,±l,.... (III.50) 

eBY 

We point out that if we imposed the periodicity condition at xq 7^ the 
allowed eigenvalues would change. Nevertheless, since at the end of the cal- 
culations we will perform the continuum limit y — 00, we are not interested 
in the exact values of the allowed eigenvalues but only in the eigenstate den- 
sity ^(xo) which is 

, di eBY 

= -^ = -7^ (III.51) 

dxo /vr 

independently of Xq. For notational simplicity, we will still indicate the 
discrete eigenvalues xq-/ as xq, then, taking into account the substitutions 
()lll.46jl - ()lll.48jl . the function ()lll.44jl becomes 

K,xoi^,y) = :^^^^)"\/^^^^P [(x - xo)^ - il/(x - 2xo)] I . 

(III.52) 

The numerical factors have been chosen in such a way that if we define 
[analogously to the two-dimensional spinors given in Eq. (IX26|) ] the 

two-dimensional spinors 



<|.S(r)^liV^/;0'„,,„(x,y) (III.53) 



exp{ikz) 
with 

J = {n,k,(T,Xo}, (III. 54) 

they result normalized as [see also Eq. ()A35|) ] 

dr^'j{r)<!>[j,{r)=6jj, (III.55) 

with 6jj> = Sn,n'Sk,k'Scr,a'Sxo,x[^. Finally, the zero order spinors uf\r) are 
given by [see Eqs. (jAljl and ()A6jl ] 

uf\r) = nlmT\^) (ni.56) 
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where 



ur\r) 



1 



(0) , rr 



\ 2w 



(0) 
R 



with 

and are normahzed as 



0s w^j^ + \/(j)t'm 
R = {n, k, a} 



dv,/^^uf^\r)U^Pir) = 5j,j. 



(111.57) 



(III.58) 



(111.59) 



Obviously, from Eq. pil.40j) the eigenvalue corresponding to this eigenfunc- 
tion is, up to first order. 



w 



(1) _ „„(o) 



Wr + 



<i)tm 



(0) 



p - 9m)— (7^ + i9E - 9p)wr 



w 



^(Ptm? + — + eB{2n + 1 + a)] <^ 1 



Xo = 

gt_ Qs fc^ + eB{2n + l + a) 
(pt (ps 0sm2 + A;2 + eB{2n + 1 + a) 
(III.60) 



where we used the definitions pi.39jl - pi.41|l of the coefficients qe, gp and qm- 
This result has a clear interpretation in classical terms. In fact, the classical 
Lagrangian of a particle in the given metric pi.4|) is 



god (x) + gll\x)v^ 



no sum 



with 



\dt J 



\dt J 



or, up to first order. 



L(^)(r, v) = -m^/4yi^^^^f)~u^ ( 1 + f 9 ^] ■ 

V 0t - (psV^ ) 



(III.61) 
(III.62) 

(III.63) 



From this expression and by defining the first-order gravitational force on 
the particle as jSHj 

f«(r,v)^^^^i^ (III.64) 
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we easily obtain that the gravitational force lies along the x-axis and its 
x-component is equal to 



Finally, if we express v'^ in terms of the momentum p defined as 
dL{r, v) 



(III.65) 



P 



mgy'{x)v 



no sum, 



(III.66) 



we have 



(pt'm? + —p^ — + 



gt , gs 



')s (pstn^ + P 



(III.67) 



where, for simplicity, we used the same symbol to indicate the force as a 
function of v and of p. By performing the obvious identification 



p2 ~ fc^ + eB{2n + 1 + a) 



(III.68) 



and by comparing Eqs. pil.67j) and pil.6fl|l . we conclude that the first-order 
correction to the eigenenergies is nothing but the gravitational potential en- 
ergy of the electron in the presence of the constant and uniform gravitational 
force (ITTLGTIi . 

In an analogous way we can write the positron eigenenergies up to first 
order as 



r,(i) 



w 



[P + eB{2n + 1 - a)] <^ 1 + 



iij^gs 



fc2 + eB{2n + 1 - a) 



/),m2 + P + eB{2n + 1 - a) 
(III.69) 

We observe that, analogously to the case of the electron eigenstates and 
eigenenergies, the quantum number nd is here substituted by xq, and it is 
pointless to keep the subscript g in the remaining quantum number Ug be- 
cause there is no more possibility of ambiguities. 

The zero-order positron eigenstates are given by [see Eqs. ()A2|) and ()A7|1 ] 



Km. 



'(0), 



with 



vrir) 



a 



2^f 




V'(r, -id) 



■XS(r) 



tin 



(III.70) 



flll.Ti: 



XS(r) 
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and [see Eq. (IX27|l ] 



X'.(r) - '-^^^^f.A,.M,y)- (111.72) 
Obviously, these states are normalized as 

j dr^^vf\r)vf{r) = 5j^r (III.73) 

and are orthogonal to the eigenstates uf\r) up to zero-order terms: 

j dr^^VP\r)Uf}{r) = j dr^^uf\r)vj?\r) = 0. (III.74) 

3.2 Computation of the one-particle states Uj{r) and 
^(r) up to first order 

We can pass now to the determination of the one-particle states Uj{t) and 
Vj(r) up to first order in the perturbation X(r, —id). From Eqs. ()III.60|) and 
pil.69jl we see that the first-order energies of the transverse ground states 
that are given by 



El + I X 

(III.75) 

are independent of the magnetic field. As we have seen in this fact gives 
the transverse ground states a particular relevance in the case of strong mag- 
netic fields because their energies are much smaller than that of the excited 
Landau levels that depend on B. In particular, we have shown that in this 
case the probability that a pair is created with both the electron and the 
positron in a transverse ground state is much larger than the other probabil- 
ities jlj. For this reason we will only calculate this probability and then we 
need to compute only the one-particle transverse ground states corrected up 
to first order. As in the previous paragraph, we will present the calculations 
of the first-order corrections to the electron transverse ground states and we 
will quote the analogous results for the positron transverse ground states. 

The first order corrections to a given zero-order electron transverse ground 
state come from the coupling of this state to the following classes of states: 

1. the zero-order electron transverse ground states with the same energy 
(all but the state to be corrected); 
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2. the zero-order electron transverse ground states with different energy; 

3. the zero-order electron states which are not transverse ground states; 

4. the zero-order positron states. 

Suppose, now, that we want to calculate the first-order corrections to the 
state labeled by Jq = {0, k, — l,Xo} = {Rq^Xq}. The states in the first class 
are labeled by {i^o, x'q} with x'q ^ xq because they have the same energy of the 
Jo-state. But, all the contributions vanish because the perturbation pil.9|l 
can not couple two eigenstates of Tt^'^^r, —id) with the same rid and two 
different Xq and x'q. Similarly, the states in the second class are characterized 
by Jq = {0, fc', — 1, Xq} with k' ^ k and, since [I{r, —id),V\\] = they do 
not give any contributions. Instead, the contributions from the states of the 
remaining two classes are, in general, different from zero then we write the 
state U^^J{r) as 

J' J' 

where the primed sum does not include the transverse ground states and 
where jSlj 

P£.' = (0) ^ (0) / dr^,U^^\r)I{r, -^^)ufJ{r), (III.77) 

QZ' = (0) ^ .(0) / dr^,V}?^\r)Iiv, -^^)ufJ{r). (III.78) 

In these equations we introduced the zero-order energies of the transverse 
ground states defined as [see Eqs. ()lll.27jl and ()lll.28jl ] 



with ^0 = {0, fc, +1}. We start by calculating the coefficients P^j^j/- From 
the expression pil.9j) of the interaction Hamiltonian and from Eqs. pil.56|l 
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and pil.57|) we have 



>(i) 

Jo, J' (0) 



-(0) 



X 



V^(^P - gM)m - + ) + ^^ef / rfr$'„t, ^, (r)a;$'o,.„(r) + 



+ 



X 



■V<h{gp - gM)m - Y + + QESk^ 



(0) ^ ^ / o?r<,^,(r)V'(r,-z5)xV'(r,-z5)<|.'o,.^,(r) 

(III.80) 



By using the orthonormal properties of the functions $^ (r) [see Eq. pil.55|) 
and Eq. ()B5|1 we obtain 



p(i) 



2e 



(0) 



X 



V2ZB 



^n',+l^o-',-l + 



-^(9P-SM)m-f«+4°V9^4' 



(0) 



X 



V2eB 



(III.81) 
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where we defined the coefficients 



d(i) 

k,xo 



ciO) _ JO) 



\ 



2S 



(0) 



X \ vMdM - 9p) 



k 

m 



A 



2e 



(0) 



V2eB 



+ 



(0) 



1 + 



(III.82) 



k,xo — 



-(0) 



X 



(III.83) 



with [see Eqs. ()lll.27j) and ^\\\.2H\f 



c(o) — (0) 



,(0) 



„r,(0) _ „r,(0) 



' 0jm2 + ^ (F + 2e5) 



(III.84) 

the zero-order energy of the ffist excited Landau level. In the same way we 
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can write the coefficients Q^j^ ji as 



a 



A(9P-<;M)m-f(4»'-4»')+9.^fl(-l, 




4>s w^2) + 



X 



X 



+ 



X 



-iSn'fiSa'-l - kXo6n'fl5a',+l 



+ 



-V^tigp - 9M)m - yl^i - w^S) + 9Eel 



(0) 




X 



k 



(III.85) 



with 



2U 



-(0) 



M9m - 9p)mJ -^kxo, 



■— I, ™_ — 



1 



4°^ + VWtm 



,(0) 



X 



X 



M9m - 9p)m - y(4°^ - ef) - gssl 



(0) 



1(1) 

p(0) , (0) 



4°' + VAm 



(III.86) 




(III.87) 



:X 



X { V<h{9M - 9p)m 



+ 



+ 



+ 



f(4"'-4°')+te4' 



(0) 



(III.88) 
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If we also define the coefficients A^^],^ as 

Ail = f^o (III.89) 

then, the first order transverse ground state U^'^l ^^^{r) can be written sim- 
ply as 

- D^') V^^^ (r) + if^^^ 1/^°) (r) - F^'^ (r) 

(III.90) 

The term ^'kJco^o%-i,xoi^) = dE^oUQ^k-i.xoi^)/'^ ^as been added to compen- 
sate for the factor (1 — QeXq) in the scalar product pi.60j) and then to have 
the states correctly normalized up to first order, as 



We observe that, even if A^.^, . . . , F^.^ are all first order-quantities in the 
couplings gE, Qp and qm, the coefficients B^^, C^^, and F^^ go, for 
strong magnetic fields, as 1/ \JleB and this circumstance makes them, in 
general, much smaller than A^.^ and D^^. We will exploit this observation 
in the next paragraph. 

Finally, with analogous calculations it can be shown that the first-order 

O,fc,+l,xo^ 



positron transverse ground state V^^,^ ^ (r) can be written as 



k,XQ *^0,fc,-l,xoV^/ 



'-^fc,xo'^0,-fc,-l,xol^/ ^ ''^k,XQ^ Q-k,^\,X(iV- ) ^k,xcy\-k,-\,XiiV- 



(111.92) 

The states y^^+ixj^^) ^^^"^ normalized as 

Kl^i;^) = ^^o,.^ (111-93) 
and they are orthogonal up to first-order to the states uf^ir): 

<^^iJ^) = K^^i^^)=o• (111.94) 

Before continuing our computation of the instantaneous eigenstates of 
the Hamiltonian H^^\t) we want to show a nice property about the mean 
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value of the velocity operator a in the first-order eigenstates pil.90|) and 
pil.92|) . In particular, we will calculate the mean value of the x-component 
of the velocity given by Vx = oix and of the component of the velocity in 
the axis orthogonal to the and to the direction of B given by vx_ = 



OyCos-i? — azsin'i? [see Eq. (IIILlll ] in the eigenstates 



and Vf\v) 

Jo 



Vr 



(1) 



0,A:,+l,a:o 



r) 



1(1- gEx)U^ 



(III.95) 



^^^^ 'v) [uy COS d — Uz sin d) u'^l _^ 

(III.96) 



and 



Jo I I Jo 

'(1) 



J rfrv^(l - ^?i.x)<^+i,.o(r)«^-<!+i,.o(r)> (ni.97) 



Jo 



mi 



- 9Ex)V^^l^^^^{r){ay cos^? - sin i3)V^ 



,(1) ( 



(III.98) 



Now, we proceed by considering only the first two mean values because that 
concerning the positrons can be calculated in an analogous way. By using the 
transformation properties ()A43jl - ()A45|l and the definitions a± = (ax^i(yy)/2 
we have 



(1) 

Jo 



(III.99) 



(1) 

Jo 



■J (Ill.lOO) 



r = 
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where, by keeping only the terms up to first order 



'(1) 

Jo 



d{1) rj(0) / N _ -(-Vi; Tj 



(1) rr(0) 



0,fc,+l,xo 



^ / k,xo l,—k,+l,xo'' 

(III.lOl) 



r + 



'(1) 

Jo 



(1) T/(0)t 



'k,xn ^1 



XQ ^ l,-fc, + l,Xo'' 
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Now, in appendix D we indicate the technique to calculate these matrix 
elements and the results are 
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then, from Eqs. ()lll.99jl and (jlll.innjl we obtain 
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and, analogously. 
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This results are just the quantum counterpart of the so-called (E x B)-drift- 
velocity effect typical of the motion of a charged particle in the presence 
of a uniform and constant electromagnetic field (E, B) [23]. In our case, 
obviously, the electric force is substituted by the gravitational force [see Eq. 
(1111.6711 ] and the sign of the drift velocity depends on the sign of the charge 
of the particle. 
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3.3 Instantaneous first-order eigenstates and eigenval- 
ues of the Hamiltonian H'^^\t) 

The results obtained in the previous paragraphs will be used here to solve 
our initial problem: the determination of the instantaneous eigenstates and 
eigenenergies of the slowly varying Hamiltonian 

H^^\t) = j dr^/^^{l- gEx)'^\T,t)n^^\r,-id,t)<iJ{r,t). (III.108) 

In fact, we have only to substitute everywhere in the results just obtained the 
constant magnetic field B with the time-dependent magnetic field B(t) given 
m Eq. (fiTT2|) . In this way, the one-particle eigenstates both at zero and first 
order will depend explicitly on time such as the eigenenergies apart from 
those of the transverse ground states that do not depend on the magnetic 
field. 

With this prescription if we expand the field ^^(r, t) with respect to the 
first-order instantaneous basis [Uj^\r,t),Vj^\r,t)] as 

nr,t) = Y,[c'^j\t)u!j'\r,t) + d'^j^\t)vi'\r,t)] (III.109) 



the instantaneous eigenstates of the Hamiltonian ()II.63|) are given by the 
Fock states 



\{nj}{ty,{njm^\c^'^\t)r' f4f(t) 



■■■\m) 

(III.llO) 

with |0(t)) the instantaneous vacuum state, while the instantaneous eigenen- 
ergies are 



(III.lll^ 



In this framework the creation of a pair at time t with the electron in the state 
labeled by J and the positron in the state labeled by J' is the transition from 
the vacuum state |0(t)) to the pair state = c^jif)d\,{t)\^{'t)) ■ 



4 Calculation of the production probabilities 

As we have said in the previous section we can calculate the probability 
of producing an e~ — e"*" pair in the presence of the slowly varying magnetic 
field flII.12jl and of the static gravitational field described by the metric tensor 
pi.4j) by means of the adiabatic perturbation theory up to first order in the 
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time-derivative of the magnetic field. In order to avoid the possible confusion 
between the "first order" relative to the adiabatic perturbation theory and 
the "first order" relative to the gravitational couplings qe, Qp and gu, in what 
follows we will always refer to the second one. In particular, the symbol 
indicates quantities that are first-order in the gravitational couplings. 

From now on, the gravitational field will not play any further role: we took 
into account its presence by correcting the one-particle electron and positron 
eigenstates and eigenenergies. Now, following the adiabatic perturbation 
theory, the operator responsible of the pair production is the time derivative 
of the second-quantized Hamiltonian pi.63jl that, by using Eq. flll.52jl . can 
be written as 



H^^\t) = / t^rv/0f(l - gEx)'^^{r, t)d^n^^\Y, -id, t) ■ B(t)^(r, t) 





gEx)^\r, t)J^[l-{gp- gE)x] | (r x a) ■ B(t)^(r, t) 



1 - 5(px)^"^(r, t) (r X a) ^(r, t). 



(IV.l) 



Starting from the second line of this equation the physical meaning of the 
result can be understood. In fact, the vector a can be interpreted as the one- 
particle relativistic operator corresponding to the velocity of the electron, 
then the quantity (r x ck) ■ B(t) = — [r x B(t)] • a is proportional to the scalar 
product of the external electric field —dA{r,t)/dt [see Eq. pi.25|) ] and of 
the electron velocity that is to the work per unit time done by the induced 
electric field itself. On this respect, we want to do a couple of observations 
that can also be referred to our previous papers |21 El E! • Firstly, the induced 
electric field E(r, t) = —dA{r,t)/dt is always perpendicular to the magnetic 
field B(t). Secondly, due to the presence of E(r, t) we may conclude that this 
last field is the responsible for the production according to the well-known 
mechanism proposed by Schwinger in |SH] • There are however two differences 
to be stressed: 

1. because of the presence of the time depending magnetic field B(t) the 
electric field E(r, t) is rotational, so it does not admit a scalar poten- 
tial and the interpretation of pair production as tunnel effect is not 
straightforward; 

2. the one-particle states of the produced electron-positron pair are states 
in a magnetic field and they very different from the one-particle states in 
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an electric field (for example, the transverse motion here is quantized 
and the energy levels with different are well separated from each 
other) . 

Now, we pointed out in the previous section that the probability that a 
pair is created with both the electron and positron in a transverse ground 
state is much larger than the others probabilities. If the electron is as- 
sumed to be created in the transverse ground state state Uj^^{r,t) with 

Jo = {0, k, —1, Xq} and the positron in the transverse ground state ^l^''(r, t) 

with Jq = {0, k', +l,a;g}, then the matrix element of the transition is given 
by 

ij«,(t) = (U(t);ij,(t)|ij«(t)|0(t)) = 




eB(t) 



j rfrv/03(i _ gpa:)U^l)\r,t) (r x a) V^{r,t) 

(IV.2) 

where \0{t)) and |1 jo(t); Ij'(t)) are the vacuum and the pair state at time t 
respectively. The creation amplitude at time t can be calculated from this 
matrix element as |2Z] 

where we used the fact that the first-order energies e^.^^^ of the transverse 
ground states do not depend on B(t) and then on time [see Eq. pil.75|) ]. 

In the following we will distinguish two different time evolutions of the 
magnetic field: 

1. rotating magnetic field: B(t) = B^(t) = B^{0, smut, cos ut); 

2. magnetic field varying only in strength: B(t) = B|(t) = {0,0, B^{t)) 
with 

B,W^B,.(B.-B,)exp(-1) ,IV.4, 

and Bi < Bf. Obviously, in both cases the magnetic field is assumed to be 
strong and slowly varying that is [see Eqs. pi.l5|l and pi.l6|l 

2 

Br. > — , (IV.5) 

e 

a; < m (IV.6) 
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and 



2 

m 



» — , (IV.7) 



e 

Bf — Ba 1 

— -- < m. IV.8 

Bi T 

The first configuration describes better that case in which the magnetic field 
is produced by a rotating body while the second one concerns a situation in 
which it is generated by a collapsing body. 



4.1 Rotating magnetic field 

In the case of rotating magnetic fields we know that in Minkowski spacetime 
it is already possible to create a pair with both the electron and the positron 
in a transverse ground state . Our task here is to calculate the corrections 
due to the presence of the gravitational field to this "fiat" creation probability. 
For this reason and reminding that for strong magnetic fields the coefficients 
A^o and D^f^l^ are much larger than B^^^^, q^l^, E^^l^ and F^^^^^, we can keep 

for simplicity only the terms proportional to A^^],^ and D^^],^ (apart from 
the zero-order term) in the expressions of the first-order transverse ground 
states.^ With this simplification the matrix element pV.2|) can be written as 

^Note that in the case of purely rotating magnetic field the coefficients A^.^],^, . . . , F^,^],^ 
are time-independent because they depend on the strength of the magnetic field. The same 
thing can be said about the electron and positron eigenenergies and about the rotated 
eigenstates U'j^\r) and V'/'"^\r). Instead, the eigenstates Uj^\r,t) and vj'^\r,t) depend 
on time through the rotation operator TZxiiSit)] — TZx{ojt). 
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The matrix element does not depend on time and, if we keep only the first 
order terms in qe, Qp and qm, we have 



H 




euBr 



/(O) 



+ / dr^ 



a(1) r/'(0)t / N 
^k,xo'^0,k,-l,xo\'- J 



r)(i) T.'(o)t / A 



(1) T/'(0) 



(IV.IO) 



Now, from Eq. ()A14|) we see that the matrix elements of the operator (1 — 
gpx)x between two transverse ground states that have both n = are equal 
to that of the operator (1 — gpXo)xo + l/{2eB) or, by coherently neglecting 
the last term, to that of {l—gpXQ)xo. In this way, the matrix element pV.lOjl 
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becomes 
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By using the explicit expressions of the transverse ground states given in 
appendix C it can be easily seen that 
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By substituting these results and the expressions (jlll.86j) of the coefficients 
1) 

;,xo 



D^^i in Eq. (llV.lljl . we obtain 
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where we used the definitions pi.39|) - pi.41|) of gE, gp and gM- Because of 
the presence of the Kronecker delta functions, the only transition amplitude 
different from zero is given by [see Eq. pV.3j) ] 



7, 



(1) 

r\;(0 ,k , — l,xo;0 ,—k ,+1 ,xq) 
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If we square the modulus of this expression and multiply it by the number 
of states [see Eq. ()lll.51jl ] 



(LNr^ = — axo X —ak 



(IV. 17) 
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we obtain the differential probability that a pair is present at time t with the 
electron (positron) between xq and Xo+dxo and with longitudinal momentum 
between k and k + dk {—k and —k — dk): 

\ (pt (ps (psm? + k'^ J 
(IV. 18) 

where the continuum limits F ^ oo and Z ^ oo have been tacitly performed 
and where the "physical" first-order volume c/y*^-'^) = (f)l^'^{l—3gsXo/(j)s){YZdxo) 
has been introduced. In conclusion, we obtain our final result by averaging 
Eq. pV.18|) with respect to time and by dividing it by the phase-space 
volume dV'^^^dk: 

/ dP^\xo,k;t) \ ^ 
\ dVdk I 27r2 04 

/ Qt Qs k"^ — 5(i>gm'^ 
V 9t (Ps (PsTn^ + k^ 

This result is to be compared to the "flat" probability that can be obtained 
by putting Qt = Qs = ^ and (pt = 4>s = ^ [see Eqs. pi.7|) and pi.8|) ] and 
that is identical to the result (1) in i7j. We note that, since in our approx- 
imations gt/(pt can be assumed to be much larger than gs/<ps, the creation 
probability increases for negative values of xq. Analogously to the sign of the 
first-order correction to the eigenenergies [see discussion below Eq. pil.6fl|l ]. 
this fact can be understood in classical terms by observing that the gravi- 
tational force pil.64|) associated to the metric pi.4|) has only the negative 
x-component pil.67j) and then, the electron and the positron are more likely 
created at < 0. Nevertheless, we point out that the first-order correction 
of the probability is not proportional to the gravitational force but this is 
due to the fact that the operator responsible of the pair creation is not the 
Hamiltonian but its time- derivative which has a completely different struc- 
ture. From this point of view, it is worth noting that the physical origin of 
the presence of the length xq in Eq. pV.lQjl is twofold. In fact, since, as 
we have said below Eq. pV.lll . the matrix element of the time-derivative of 
the Hamiltonian is proportional to the matrix element of the power supplied 
by the external electric field which grows linearly with the transverse spatial 
coordinates, then the probability ()IV.19|) results proportional to Xq. Instead, 
the first order correction due to the gravitational field is proportional to Xq 
because the gravitational force is uniform [see Eq. pil.65j) ]. In this sense. 
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the dependence of our final results on the quantity Xq is unavoidable since it 
is a direct consequence of the (approximated) model we used that is of the 
assumed uniformity both of the magnetic and of the gravitational field at a 
microscopic scale. A more realistic model should take into account that these 
fields are not actually uniform in the whole space and that they vanish far 
from their source but the calculations would be too complicated. Concerning 
this point, we want to stress that from the derivation itself of Eq. pV.19|) . 
a clear physical meaning can be assigned to the quantity xq. In fact, it can 
be interpreted as a typical length scale in which the magnetic (gravitational) 
field produced by the astrophysical compact object can be assumed to be uni- 
form and this fact gives the possibility to do also a quantitative prediction 
from Eq. (I1V.19I1 . 

A final note concerns the presence in Eq. pV.19|) of an overall factor 
4>s/4>t > 1 that increases the creation probability with respect to its "fiat" 
value. 

4.2 Magnetic field varying only in strength 

We have seen in that, if the magnetic field varies only in strength it is 
impossible in the Minkowski spacetime to create a pair with both the electron 
and the positron in a transverse ground state. Instead, we want to show here 
that this process is allowed in the spacetime metric pi.4|) because of the 
corrections induced by the gravitational field on the transverse ground states 
of the electron and of the positron. In fact, if we choose the same electron 
and positron states used in the previous paragraph, the transition matrix 
element (IIV.2jl in this case becomes 

(IV.20) 

where we used the fact that, since the magnetic field lies in the z-direction 
for every t then TZxl'^it)] = I. 

Now, the selection rule (50) in p3 concerning the spin of the transverse 
ground states allows us to conclude that for the zero-order transverse ground 
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states the following equalities hold 



drU'jf\r,t) (rxa)^\/f)(r,t) = 0, 



(IV.21) 



J drUjf\r,t) (r x cx)^Ujfir,t) = 0, (IV.22) 
/ rfrVf (r, t) (r X a)^ ^ (r, t) = (IV.23) 

J -^0 Jq 

where Jq = {0, /c', — 1, Xg} and Jq = {0, fc', +1, Xg}. By exploiting these 
equations and by keeping only the terms up to first order we can write Eq. 
(jrvTWijl as 
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This expression can be further simplified if we use the definitions ()A14|) and 
()A16j) of the operators x and y. In fact, by using the substitutions pil.46|l - 
pil.47jl the operator (r x ot)^ can be written as 



(r X a)^ = xoy - ya^ 
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V2e5 

i 



{a + a"^) 



[a — a 



i{a- — «+) — 
q;_ + a+) - 



Xq + iyo + \ —^a 
eB 



ia+ I Xo - iyo + \ j —a ] . 



(IV.25) 

Now, from the expressions ()C1|) and ()C2|) of the transverse ground states we 
conclude that only the operator ia;_(xo + iyo) gives a non- vanishing contri- 
bution in the first integral of Eq. pV.24j) and that, analogously, only the 
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operator — ^^^(xo — iyo) gives a non- vanishing contribution in the second 
one, then the matrix element ()IV.24|) can be written as 




where 
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and where we used Eq. (jA23|l . The calculation of these matrix elements 


and pV.28|l is quite tedious and the technique to perform it is sketched 
in appendix D. In particular, it can be shown that 

(i'rwi-it^rw)=-(<'wi".it^rw)= 



(t)t{9M-9p)mk6k-k'S X 



(IV.29) 

and, for this reason the terms with the derivative with respect to Xq cancel 
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each other and the final expression of H^^^^^j^j^^{t) is 
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We want to point out that the disappearance of the electron charge — e is only 
a consequence of the strong field approximation B^{t) ^ B^r- In particular, 
it would be incoherent to conclude that this matrix element does not vanish 
in the limit e ^ because in this limit Bcr = m? /e — * oo and the strong 
field condition can not be satisfied. Nevertheless, we will see that the creation 
probability which is physically meaningful will contain e. 

The creation amplitudes at time t can be calculated by means of Eq. 
pV.3jl and the only one different from zero is equal to 
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where only the first-order terms in the gravitational couplings have been 
kept. Now, by substituting the expression pV.4j) of the magnetic field B^{t) 
and by calculating its time derivative 
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the creation amplitude can be written as 
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with s' = t'/r a. dimensionless variable. Now, in our model r is a macroscopic 
time connected to the typical evolution times of a neutron star or of a black 
hole, then we can safely assume that e^^^r ^ 1. This allows us to give an 
asymptotic estimate of the remaining integral in the limit t ^ oo. The result 
is 
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Finally, by squaring the modulus of this expression and by multiplying it by 
the number of states at t oo 

' 27r " 271 (27r)2 ^ ' 

with = ^^/^Y Zdxo the "physical" quantization volume up to zero order 
[see Eq. pi.5|) ]. we obtain the differential probability that a pair is created 
with the electron (positron) between xq and Xq + dxo and with longitudinal 
momentum between k and k + dk {—k and —k — dk) 

(1V.36) 

In this equation the continuum limits y ^ oo and Z ^ oo have been 
performed and, since the probability is already proportional to Qs and our 
calculations are exact up to first order in Qs and gt-, it is enough to use the 
zero-order "physical" volume dV^^\ Finally, the corresponding probability 
per unit volume and unit longitudinal momentum is given by 



dP\'^\xo,k]t (X)) f eBf f Bf - Bi 



^mxok 



dV(°)dk \l6Ti(psJ 0t (0.m2 + A;2)3 \^ rfi,- 

(IV.37) 

As in the case of a rotating magnetic field we observe the presence of the 
overall factors 4>s/4>t and x^. Also, in [2] it was calculated the total probability 
that a pair is created in a slowly varying magnetic field with fixed direction 
but in the Minkowski spacetime [see Eq. (20) in this paper]. We remind 
that in this case either the electron or the positron must be created in a 
state which is not a transverse ground state. In order to have a quantity to 
be compared to the total probability given in Eq. (20) in [2], we have to 
integrate Eq. ()IV.36|) with respect to Xq and to k. After these integrations 
it can easily be seen that the order of magnitude of the ratio between our 
probability and that given in [2^ is^ 



1 f9, 
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If we substitute the expressions pi.7|) and pi.8|) of (pt, 4>s and gg we see that 
because of the inequality pi.lO|) and of the strong- field condition eBf ^ 

^We have to be satisfied of an order-of-magnitude comparison because the time evolu- 
tion of the magnetic field used in is different from that used here. 
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then rj results much less than one and then the gravitational effect is in our 
approximations small. Nevertheless, the effect is there and it is reasonable 
to imagine that it can be amplified in the presence of a real gravitational 
field which is not restricted by our assumptions. For this reason we decided 
to consider also a particular case slightly different from the one at hand and 
where the gravitational field can be treated without approximations and this 
is the subject of the next paper. However, it must be pointed out that the 
perturbative approach gives the possibility to treat more general magnetic 
and gravitational field configurations while the strong field case the two fields 
must be necessarily parallel. 

5 Conclusions 

In this paper we have modified our previous results given in [2 El El about 
the production of — pairs in the presence of strong, uniform and slowly 
varying magnetic fields taking into account the presence of the gravitational 
field represented by the metric tensor pi.4|) . We point out that this metric 
tensor resulted from some approximations that we have done in section (j2)) 
and that we resume for the sake of clarity: we have neglected the gravitational 
field due to the magnetic field energy, we have neglected the gravitational 
effects of the angular momentum of the stellar object and we have assumed 
to be not to close to the event horizon of the body. 

We have treated the gravitational field perturbatively up to first-order in 
the couplings qe, gp and qm and we have seen how its presence modifies the 
one-particle eigenstates and eigenenergies of the electron and of the positron. 
This circumstance, obviously, changes the pair production probabilities we 
calculated in absence of the gravitational field, i. e. in the Minkowski space- 
time. 

In particular, we have reexamined the case of a purely rotating magnetic 
field jl] and of a magnetic field varying only in strength |2]. Firstly, in both 
cases we have found that the production probabilities contain an amplifica- 
tion factor > 1 [see Eqs. ()IV.19|) and pV.37|) ]. Also, in the case of the 
purely rotating magnetic field we have seen that the presence of the gravita- 
tional field in the plane orthogonal to the magnetic field brakes the rotational 
symmetry in this plane and makes more likely that a pair is created at xq < 
see Eq. fllV.1911 ]. Instead, in the case of a magnetic field with fixed direction 
we obtained a qualitative new result: in the presence of the gravitational 
field it is possible to create a pair with both the electron and the positron 
in a transverse ground state. Nevertheless, we have seen at the end of the 
previous section that this probabihty is a small quantity with respect to the 



44 



total probability that a pair is created in the Minkowski spacetime. But, 
this result is a consequence of the fact that the gravitational field has been 
treated perturbatively and this treatment could not fit what it happens, for 
example, near the event horizon of a black hole. From this point of view, 
the information we have gained is that in the presence of a gravitational field 
this new effect is there. In the next paper we analyze a particular case in 
which the gravitational field can be treated without approximations and we 
will show that in that case the effect is relevant. 
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Appendix A 

In this appendix we want to quote some properties of the eigenstates Uj{r) 
and Vj{r) of the zero-order Hamiltonian pil.8|) . Firstly, we remind that 
j = {nd,k,a,ng} embodies all the quantum numbers and that Uj{r) and 
Vj{r) can be written in the form |3] 

«,(r)=7^t(^^(r), (Al) 

v,{r)=nl{^mr) (A2) 

where 

7^,(^9) = exp{-i^J^) (A3) 

with J7x = Cx + ^x the x-component of the one-particle total angular momen- 
tum operator and t? is defined in Eq. flIII.3|) . The rotated spinors u'j{r) and 
v'Ar) are the solutions of Eqs. ()III.21|) and pil.22|) with the vector potential 



and 



A'(r) = -i(rxBO (A4) 



B' = I 1 = 1 I 1 (A5) 
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[seeEqs. (ITTLT11 - (ITTL4|1 



, and are given by |1] 
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.7,(0) 




0s wf^ + y/cfhTTl 

X'{r) 



(A6) 



(A7) 



where the indices r and q and the energies Wr""* and Wq^^ have been defined 
in Eqs. fllll.24j) and (I11L25D and in Eqs. (1111. 27jl and (llll.28j) respectively 
and where the numerical factor 1 / -y/0f has been inserted to make the spinors 
correctly normalized with respect to the scalar product pi.60jl at zero order: 

j dv^^u]ir)uAr) = j dr^,u';{T)u'j,{r) = 5,,,, (A8) 
I dr^,v]{r)vj{r) = j dv ^^v'^ {v)v]{v) = 5,,, (A9) 



with = 5,, 



Now, we give a more explicit expression 



of the operator V'{r,—id) and of the two-dimensional spinors and 
Xj(r) appearing in Eqs. ()A6|) and ()A7|) . To do this, we need to express 
the operators x, —id/dx, y and —id/dy with respect to the destruction 
and creation operators (that we indicate as ad, aj^, ag and respectively) 
corresponding to the quantum numbers Ud and Ug and vice versa. These 
expressions can be found in [BBI and are given by 



ad 




2 

eB 




2 

eB 




2 

eB 




{x - iy) + 
(x + iy) - 
{x + iy) + 
{x - iy) - 





2_ /d_ 

eB \ dx 

2_ f d_ 

eB \ dx 

2_ f d_ 

eB \ dx 

1_ 

eB V dx 




ld_ 

i dy 

ld_ 

i dy 

ld_ 

i dy 

ld_ 

i dy 



(AlO) 
(All) 
(A12) 
(A13) 
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and by 



1 / 2 , t U 1 



1 d _ 1 eB 

i dx 2i 

1 fY 

1 [eB 



a\ + ad- a\) 



eB 



ad + a\), 
1 /e5' 




y 



(S - 4 - ttrf + ai) 



1 



yo 



iV2eB 



7^ = -2V^^"^ + "^~"^~"^^ 

where we introduced the operators 



eB 1 



2 




+ a 



(AM) 
(A15) 
(A16) 
(A17) 



(A18) 
(A19) 



Starting from the well-known commutation relations among the operators x, 
—id/dx, y and —id/dy it is easy to see that 



1, 





and that 



'eB' 



[a;o,l/o] 

In particular, this last commutator allows us to put 

1 1 d 

yo 



eB i dxo 



(A20) 
(A21) 

(A22) 
(A23) 



Now, the operator V'(r, —id) in Eqs. ()A6j) and ()A7jl is defined as [see Eq. 
(21) in I?] 

V'(r, -td) = a- ■ [-td + eA'(r)] (A24) 
then by substituting the expressions ()A14jl - ()A17|l it can easily be shown that 



V'(r, -id) = i^2eB{a\a- - ad<J+) + <^z--g^ 



(A25) 



where a± = {ax ± iay)/2. It is also convenient here to express the two- 
dimensional spinors <^'j{r) and Xj(r) in Cartesian coordinates instead of in 
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cylindrical coordinates as in Eq. (22) of Ref. |3]: 



(A26) 
(A27) 



In these expressions Z is the length of the quantization volume in the z- 
direction, 



exp(-iA;2;) 

^j^"^^ = — — f-^^n„nAx,y)- 



fU 



f'-l 



and where the scalar functions 



eB 1 1 



27r h\W. 



ai)'Ha;)'^ exp 



eB{x^ + y 



2M 



(A28) 



(A29) 



depend only on the transverse coordinates. In Eqs. ()A26j) and ()A27jl the 
operators a\ and are supposed to be expressed as in Eqs. (jAllj) and 
()A13|) and from those equations it can be shown that 



(ap"» exp 



(a^)'''' exp 



2M 



eB 



{x — iy)""" exp 



eB{x'^ + y 



2M 



(A30) 



eBY" 

— j (x + iy)"'' exp 



eB{x^ + y^ 



(A31) 



and then that the transverse functions ()A29|) can be written as 

' ' eB{x^ + y^) 



eB_ 



[x — iy) ^ exp 



(A32) 

Finally, we want to derive some orthonormalization properties of the two- 
dimensional spinors 9?j(r) an d x ^ ( r) and of the spinors Mj(r) and Wj(r). From 
the commutation relations ()A20|) it can easily be shown that, given two 
quantum numbers and n'^ with n'^ > n^, then 



rid 
.1=0 



(A33)' 
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and, symmetrically, 



ia,r'(al)"'. = (a AT' (air'.-''' = (l+ala,)"' (a\f.-"' 



.«=0 



nr. . , , 



(A34) 

if > Ug (if < Ud or < rig analogous relations are obtained). With 
the help of these two equations and by reminding that [a^, a^] = it can be 
seen that the following orthonormalization relations hold: 



J'' 



(A35) 
(A36) 



Now, if we calculate the square of the operator V'(r, —id) as given in Eq. 
(I A 2511 we have 



+ 2eB 



+ eB {2a\ad + l + a^ 



(A37) 



then, from the expressions ()III.27|) and ()III.28|) of the zero-order eigenener- 
gies, we obtain 



^^V^/W ..(0) , /X- (0) , f-r ^i'(^) - (0) 



Wr + y/ (f)tm Wl./ + ^/ (j)tm 



X.,'(r) = -^j^ — ^^3,3'- (A39) 



Finally, from these relations and from Eqs. ()A1|) . ()A2|) . ()A6|) and ()A7|) we 
derive immediately the equalities 



(0) "JJ'' 



(A40) 



where we used the fact that [see Eq. flA3|l ] 

7^,(t9)/37^t(^9)=/? 



(A42) 
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By the way, since they are often used in the main text, we also remind the 
following transformation properties of the matrices a 

n,{^)a,nl{^) = a,, (A43) 

n^^{^)aynl{^) = tty cosi? + a^sm'&, (A44) 

7^^(^?)Q;^7^J;(^?) = -ay sin '& + a^cosi? (A45) 

and of the position operator r 

n^{^)xnl{^) = X, (A46) 

n^{^)ynl{^) = ycos^ + zsin^, (A47) 

7^^(^?)^7^i(^?) = -y sin + ^ cos t?. (A48) 
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Appendix B 

We want to show explicitly here that 

1r.,n,;r+,n'^= j dv ^/^^ul_ ^^^{r)I{r , -ld)Ur+,n'^{r) = , (Bl) 

^r+,n,;r.,n'^= j dr ^/^^ul^^^^{r)I {r , -id)ur_y^{r) = (B2) 

where r_ = {ud + 1, k, —1} and r+ = {ud, k, +1}. We will prove only the 
first of these equalities being the other analogous. From Eq. pil.9jl we have 



(0) 



dr^^u'^_ (r) ^/(jhigp - gM)l3m - {gp - gE)w. 



(B3) 

where we used Eqs. flA43|) and ()A46|) and the fact that w'^} = wi^^ . This 
equation can be rewritten in the form [see Eq. ()A6|) ] 



X, 



(0) I /jt 
Wr_ + ycptm 



r-,ng;r+,n- 



2w 



(0) 



^Mgp - gM)m - {gp - gE)w^r} / dYip'}_ {Y)xip'^ 



+ 



^t{gp - gM)m - (gp - gE)w. 



(0) 



X 



— / dYip;_ 



(0) , ^ (0) , ^^ 'fr^^n'S^) 



Wr_ + ^/4>t^m Wr_ + V0tm 

(B4) 

Now, by observing that the spin of the two states is different we realize 
that the first integral vanishes because of the orthonormalization relation 
()A35|) . Also, from the definition ()A24|) of the operator V'(r, —id) and from 
its expression ()A25|) we have that 

V'(r, -id)xV{Y, -id) = [V'(r, -id), s]V'(r, -id) + xV^{r, -id) = 

= -ia^V'ir, -id) + xV'2(r, -id) = 



(B5) 



+ X 



+ eB {2a\ad + 1 + a, 



where we used Eq. (lA37j) . The only term of this operator that can change 
the spin of the electron is the second one, but it can not change the value of 
the quantum number n^. This means that the second integral in Eq. ()B4|) 
also vanishes and this completes the prove of Eq. ()B1|) . 
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Appendix C 

We want to give here the exphcit expression of the zero-order transverse 
ground states and of the zero-order states corresponding to the first-excited 
Landau levels both of the electron and of the positron. We remind that 
we used them to calculate the transition matrix elements pV.ll|) . piI.101|) . 
(I111.102j) . (II V. 271) and (llV.28jl . These states can be easily obtained by substi- 
tuting Eq. |A25|) in Eqs. ()lll.57j) and (jlll.Tlj) . We give only the final results 



52 



for the rotated states because, actually, we used only them 



r/'(°^ (r) 



A 



( 



2e 



(0) 





1 / 



exp{ikz) 



(r) 



f/'^°^ (r) 



r/'(°^ (r) 



(r) 

^O,k,-l,xo\'- ) 



'(0) f 
k,+l,xo ^ 



1 



^k^ + V<hm 



2e 



(0) 



\ 
I 



\ 24 



(0) 



(0) 



A 



+ v^m 



2E 



(0) 



4°^ + 



V 

/ 

V 



(0) 











0o,xo(a;,2/) 




(CI) 

\ 

exp(— z/cz) 



I 

(C2) 



exp(iA;2;) 



(C3) 

\ 



0'i,.o(a;,y) 



exp(zA;z) 



-%,xoi^,y) 



(C4) 



/ 

(C5) 

iV2^e'o,xoi^,y) ^ \ 

<^«4°'+v^"^V ~ke[x^{x,y) 


^'i,xoi^,y) ) 

(C6) 



1 

/ 




exp( 


-ikz) 


V 




exp( 


-ikz) 





In these expressions we have used the definitions pil.53j) and pil.72j) for 
the two-dimensional spinors $j(r) and Xj(r) and the definitions pil.79jl and 
pil.84|) for the energies ef^ and E^^^ . 
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Appendix D 

In this appendix we will give an example which explains how to calculate the 
matrix elements (jlll.inijl . ()lll.in2jl and ()IV.27jl and ()IV.28jl . In particular, 
we will calculate the integral 



,fc, — l.XQ * 



(Dl) 



We first observe that, given a generic spinor 5", 



5* 



S2 
S3 



(D2) 



then 



a-S 



f 














1 




















/51 \ 

^2 
^3 
V^4/ 



/ \ 

^3 



\s, j 



(D3) 



In this way, from the expressions ()(^3j) - ()(^6|l of the zero-order electron and 
positron eigenstates we see that 



a. 



I 



2£ 



(0) 







(1) 



(0). 



4itm 



Ed) 



r(l) 



-fcE' 



(1) 



/2eiJF' 



(1) 



7(0)_ 



1/^ 
exp(iA;z) 



0fm 



(D4) 
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By using the expression I^Cl^ for Uq'^J ^^{r) we obtain 



2e 



(0) 



A 



4°^ + v^m, 



2^, 



(0) 



■X 



X 




k,xo 




k 



(0) 



0tm 



(1) 



(-(1) 




t '^^k,xo 



(0) 



0tm 



(D5) 



At this point we have to substitute the coefficients B^^^^, C^^^^, E^.^^^ and 

V^^l^ by means of Eqs. (I111.82|) . (1111.83^ . (I111.87jl and (llll.88jl and many 
terms cancel each other. After some algebra we obtain 



-(0) 



(t)tm£^ 



(0) 
k 



2e 



(0) 



2S, 



(0) 



■X 



X 



X 



+ 



X 



(0) ^(0) 




V^(^M - 9p)m + — {Sf^ + ef) - gEsl 



1 



(0) 



0(0) , _(o) 




1 



(0) 



(j)tm 



c-(o) _ ^(0) 



X 




k 



c(0) , _(0) 

"T" ^k 



X 




(D6) 
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Finally, if we use the equalities 



we have 



1 + 



4°^ + V(Ftm' 



2e 



(0) 



(D7) 
(D8) 




0s 24°)4°) 



X 



cio) _ JO) 



X 



X 



(0) 



+ 



X 



ciO) , ^(0) 
^k ^k 



\/^tm{gM - 9p) \f^tm - y - ) + 9Eei 



(0)' 



2e5 




{9m - 9p)^j^ + (^P - 9E)ef 

^k 



m"^ + k"^ f 9i ^ 9s 



2eB 



k' 



m? + k"^ 



(D9) 



that is one of the equalities in Eq. ()lll.in8;i . We say again that all the other 
matrix elements can be calculated with the same technique. 
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